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ABSTRACT 



In Part I the primary longitudinal stability derivatives of a 
"conventional" airplane and a "typical" large, high speed, single 
rotor helicopter are simply compared, and fundamental similarities 
and differences discussed. Physical interpretation and meaning is 
given to various aspects of the character of the two types of air- 
craft. 

In Part II the influence of rotor flapping (tip path plane 
orientation and coning) on the stability derivatives of the whole 
helicopter is investigated, and then various stabilization schemes 
in which the flapping of the rotor is the basis for the feedback 
signal are explored. In particular are the effects of Oehmichen 
and Delta Three pitch-flap couplings on the static stability of the 
helicopter. It is found that of the techniques investigated and 
within the limits of validity of this work, only negative Delta 
Three (positive pitch-flap coupling) is effective in improving the 
longitudinal gust response characteristics in high speed flight. 
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INTRODUCTION 



This thesis is in two parts. The first is a rudimentary comparison 
of the longitudinal gust response characteristics of a "conventional" air- 
plane with a "typical" large, high speed, single rotor helicopter. The 
second part consists of the investigation of various forms of helicopter 
stabilization techniques, especially those using rotor flapping as the 
sensed element. These Wo subjects, while related, are themselves topics 
of considerable reach, and their being treated together in a master's 
degree thesis is indicative of the depth they are explored. There is a 
reasoning, however, behind this overview type of approach. 

Helicopter pilots have always complained about the so called flying 
qualities of their craft in gusty weather. "They blew around too much." 

In order to find out why they do so, and in what way, it is expedient to 
relate the problem to one whose characteristics are better known - the 
airplane. The problem is commonly formulated as a set of differential 
equations with constant coefficients, and the characteristics of the 
problem are the signs and magnitudes of the coefficients - the stability 
derivatives. Trends, order of magnitudes, and signs of these derivatives 
indicate fundamental differences and similarities between the Wo types 
of aircraft and lead to insight as to why the helicopter responds more 
severely to gusts. The rudimentary form of the analysis is such that 
only these types of conclusion can, and are, drawn. 

The first part is thus a review of the derivation of (what are con- 
sidered to be) the important stability derivatives for aircraft gust 
response, and the reasoning behind several approximations to the full 
equations of motion with the author's Wo cents worth of physical 
interpretation and understanding added. It is, so to speak, the state- 
ment of the problem to be solved in the second part. 
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In the second part, again, there is a good deal of re vie-;, the 
emphasis being on the physical interpretation and consequences of the 
findings. The influence of rotor flapping on the stability derivatives 
of the whole helicopter, and hence the equations that govern how the 
pilot will be bounced around, is investigated. The impetus for this is 
to try to discover, in a simple, fundamental manner, whether or not 
stabilization techniques involving the reduction of rotor flapping 
as a means to improve helicopter flying qualities are as "good" as they 
might seem from a more cursory glance at the problem. This is done by 
first deleting the effects of flapping in the expressions for the 
stability derivatives, and then by investigating feedback techniques 
such as Oehmichen and Delta Three pitch-flap couplings. Much time is 
spent with these two schemes because, as is shown, they form the 
basis for several other stabilizer techniques which have been proposed. 
Finally, all of these results are contrasted with a conventional body 
mounted rate and attitude SAS (stability augmentation system). 

All of the analytical analyses are complimented by dynamic analog 
computer responses. The data which is used to generate the helicopter 
models is presented in Appendix I, and serves to characterize the large, 
high speed, single rotor helicopter. Additionally, in the appendices are 
given a treatment of axis system conversion, and the definition of gust 
alleviation as it is used in this work. 
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DESCRIPTION OF ANALYSIS 



The equations of motion used as a model of the aircraft in this 
work are the standard, small perturbation equations commonly used in 
stability analyses. They are extremely simplified by the approximations 
and assumptions that form the basis for this study - and result in a form 
of equations in which only general trends and important characteristics 
can be discovered. The justification for all of the simplification is 
an interest in finding order of magnitude trends that are characteristic 
of fixed or rotating wing aircraft, rather than the details that may 
vary considerably among aircraft of the same type. Some of the assumptions 
are purely for simplicity and others are assumptions that have been shown 
to be valid for an analysis of this sort, provided certain restrictions 
are kept in mind. In particular are the effects of high gain feedback 
designs involving the rotor in the feedback loop. In such cases known 
stability boundaries are avoided. Throughout the analysis the simplifying 
assumptions almost invariably are optimistic with regards to stability. 

For this reason it is considered desirable to obtain results of a simple 
analysis before a very detailed study is conducted. 

The approximations, in addition to the standard ones (ref 10^ are: 

1) Forces, moments, and blade motions are quasi-static linear functions 

of the three aerodynamic variables A.w,^e , andMj, and the two control 
inputsA © 0 and . Implicit in this assumption is: 

a) the neglect of unsteady aerodynamics, or aerodynamic lags, such as 
the familiar down wash lag, and 

b) that rotor plane dynamics can be described by additional, redundant, 
equations in the helicopter variables. 

2) There is considered to be a plane of longitudinal symmetry, and the 
motions in this plane are uncoupled from the lateral-directioi al motions. 

3) The gusts are modeled as moving air mass fronts along body axes. 

4) No attempt is made to simulate gust alleviation, gradual gust 
penetration, or gust spatial distribution. Gust alleviation , as it is 
used in this work, is defined in Appendix 3V. 



The equations are derived according to the procedure in Reference (ll'i. 
The aerodynamic forces in the form of Taylor series, expanded in the three 
aerodynamic variables along the X, Z, e,nd about the Y axis are summed with 
the control and inertial forces and moments in those directions and divided 
by the mass and the moment of inertia respectively. The terms are 
linearized for small perturbations and the sum set equal to the trim, or 
initial vaJ-ues. This leaves equations of the form: 



XlA AlA Xvj A^a) ♦* "V “ Xb, "V )( k 0 * 

( 1 ) ~ ~ 2 o\ &{A — Xw AuJ “ "ZTo AQ — *2 ts j | + “£ & * o 

-rAvaooJ ® Ha, 6B, -v 



The stability derivatives are dimensional derivatives with the units of 
l/sec for the force derivatives and l/sec-ft for the moment derivatives. 
They are defined as follows : 

y - J__ b X - J_ XV 1^0 ~ t" \t\ 

where the partial derivatives indicate that four of the five variables 
are held constant while the derivative is taken respect to the fifth. 
This does not imply, however, that any dependent variable such as blade 
flapping, as it is treated in this work, is held constant during the 
process . 



The equations describing the motion of the rotor plane (blade 
flapping) are, consistant with the assumptions made, independent 
redundant equations of the form: 
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This work pertains to the gust response of aircraft and so a 
suitable way of describing the gusts, especially for use on the aralog 
computer and for ease of physical understanding, is used. The aerodynamic 
variable, &w, is used to represent aircraft vertical velocity with respect 
to the local air mass such that 



■=- &oc 



angle of attack 



A gust, wg is modeled as the velocity of the local air mass with respect 
to the stationary atmosphere, and the aircraft velocity with respect to 
the stationary atmosphere is T *^a. The sum is the aircraft velocity with 
respect to the local air mass such that w = (awg + Awa). 

This substitution is made into equation (l), noting that the w term in 
the Z force equation is an inertial acceleration of a body, not an 
acceleration of the air mass, s.nd should really be w^ . Ail of this is 
shown in the figure which follows . 




The equations can now be written in a form convenient for putting on 
the analog computer because the gusts are considered as separate inputs, 
or as forcing functions. They can be arranged as follows ml patched into 
the computer in a standard manner. 



G 



- ^U-CX 



o 



XvM 






-Zw 





The left hand side represents the accelerations on the aircraft due to 
the right hand side. The analog computer diagram is presented as 
figure (l). The variables printed on the diagram are really "computer 
variables" because they are normalized to maximum values. The analog 
computer responses are labelled in dimensional units and all, except 
where noted, were made with a 7y ft/sec gust input. The values for 
all of the potentiometers were derived from the data from which the 
derivatives of Appendix I were taken. 

It should be noted here that an analysis of this form differs 
from that where the transient response of the aircraft due to a control 
pulse input is studied. A control pulse input will cause different 
relative accelerations, along and about the three axes, from a gust 
disturbance, resulting in different aircraft motions. This can be 
seen as the difference between the numerators of the transfer functions 
for motions resulting from a collective input as opposed to a vertical 
gust input. 




AO 



SjAvj 



A0o 



A00 




Y 



The character of the motions are the same, but the actual responses are 
dependent on the numerators, When dealing with flying qualities the 
actual response is often of more importance than the character of the 
motion. 

Much of this analysis will consist of the interpretation of 
analog computer runs, in contrast to the solution of transfer functions 
or mode ratios and mode shapes. The a.nalog computer approach is taken 
in order to get a grasp on the basic nature of the responses and because 
of the ease in the analysis. As such it will be physically intuitive to 
keep in mind the simplified computer diagram which follows. 
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In this way one can see directly all of the parameters making up the 
airplane accelerations and, of importance in interpreting the modes 
of motion, can see the number of integrations in each loop. It is 
convenient to think in terms of electrical signals, adding arid sub- 
tracting from one another and being integrated. In the second part 
of the thesis various feedbacks are investigated. Each is merely 
a combination of the helicopter variables summed together and called 
the feedback signal. This is shown in the full analog computer di- 
agram of figure (l) where all of the different feedbacks are shown 
on the right side of the figure. Selection of one or more of the 
feedbacks is made by throwing the proper switch , which sums together 
a particular combination of all the variables and amplifies, attenuates, 
or integrates it according to the method chosen. The signal is 
then put into either the B, s ^ j p or Go Ajp amplifiers on the left 
of the diagram as a control input. 

Approximations to the equations of motion will be used frequently 
for ease of analysis. It will be shown in the first part that the 
approximations used are valid for the purposes of this work. This will 
be shown by comparing results of successively broader approximations, 
and, at the same time, will bring out the major differences between 
helicopter and airplane gust response characteristics. 

It is necessary to have a standard of comparison for a study such 
as this. Unfortunately, there is no known single parameter that en- 
compasses all of the important flying qualities parameters. Of even 
more concern is the lack of knowledge as to what the important flying 
qualities are. For the purposes here, however, rotational and normal 
acceleration magnitudes, as well as dynamic stability characteristics, 
are used as a general index of comparison. The acceleration magnitudes 
are important for structural and ride quality considerations, and the 
dynamic stability is important for control and flying qualities. 
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PART I 



ANALYTICAL AND QUALITATIVE COMPARISON OF AIRPLANE A?ID HELICOP. j XoT 
RESPONSE CHARACTERISTICS. 



In this part different approximations to the full three degrees of 
freedom of motion of an aircraft are considered. The value of each 
approximation is shown and the important points it reveals are discussed. 
As each approximation introduces stability derivatives into the equations 
they will be simply derived for both the helicopter and the airplane. 
Typical values for the derivatives are taken from data presented in the 
appendices , and in the references. 



In order to gain an intuitive understanding of the various factors 
affecting an aircraft T s response to a vertical gust, it is helpful to 
start as simply as possible. Consider an aircraft, or any body with seme 
lift curve slope represented by the stability derivative Zw, restrained 
from moving in any direction - bolted do wn in a wind tunnel for example. 
If the body is exposed to a vertical component of air velocity, Wg, then 
there is an aerodynamic force created on the body in the Z direction, and 
it .could be measured by the wind tunnel balance. The force is 

- *2 

& W 

Divided by the weight gives the normal load factor, 





A.U5 






If high frequency and, in the sense used here, higher order dynamics 
are neglected (e.g. blade dynamics, structural response), then this is the 
initial vertical gust response of an aerodynamic body. The higher order 
dynamic response is referred to as gust alleviation and effectively works 
so as to reduce the magnitude of the initial response, fhe load factor 
transfer function may be thought of as 
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where G(s) is a dynamic transfer function which, for very low frequency 
inputs can be approximated by unity. (Appendix IV has a further discussion 
of gust alleviation.) This approximation then becomes a static approximation. 
It seems obvious that the derivative, Zw, will be important in the gust 
response characteristics of aircraft. 



For a conventional airplane it can be shown (reference 11) that 



Zw 



Uc Cl 



Uo 0. 

«£/<• c 



airplane 



where is airplane relative density and a w is the slope 

of the lift curve for the whole airplane. Zw is proportional to airspeed 
but inversely proportional to size and relative density. Small, light, 
high aspect ratio airplanes flying at low altitude and high speed have 
large Zw's (order of-2) and heavy, high wing loading fighters flying 
slowly in approach have low Zw ' s (order of -0.2). 



For the helicopter, Zw is markedly different. If all of the lift is 
considered to be from the thrust of the rotor, 

W ^ ’fcvW ' AA ( C>UJa/jP 

In terms of the rotor parameters Zw is proportional to 



Q.R h CT _ 

UJ 




&y _ 



[• 



aer 

M 
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The induced velocity parameters 




are taken from reference (12) where they are shown to approach unity, 
zero, and zero, respectively, as forward speed increases and induced 
effects become negligible. Using these high speed approximations and 
noting tha,t, in level flight the second term will be negligible with 
respect to the first, then 



which is similar in form to that of the airplane, but is independent of 
airspeed. For a helicopter Zw approaches a constant value, for a given 
blade loading and RPM, at some high enough forward speed that induced 
effects are not important. It can be seen in the data of Appendix I 
that the speed range considered herein is not yet high enough to warrant 
this approximation because Zw for each helicopter is still increasirg 
even at 180 kts. It is also interesting to note that the helicopter 
represented by the circle and the square are the same helicopter where 
the larger Z w corresponds to an empty weight condition and the lower Zw 
to a heavily loaded condition. The low "triangle" values are for another 
similar helicopter which is at its maximum gross weight. 



cK. r 



acr 



for high speed. 






Dimensionalizing yields 




helicopter (high speed) 



1? 



In order to contrast airplane and helicopter Zw's it must be noted 
that Zw is highly dependent on the flight condition for a given airplane , 
whereas for the helicopter Zw approaches a constant with speed. T he 
design flight condition thus enters the picture when a comparison is made. 

It is clear that at some speed and altitude (dynamic pressure') the heli- 
copter and airplane will have the same Zw. Figure (2) indicates that, 
for six dissimilar aircraft for which data was readily available (ref. 10'), 
the speed where this happens may be near the design condition for the 
particular airplane. The values of Zw for all six aircraft near their 
design flight condition were within the range of data for four similar 
helicopters at high forward speed. For purposes of comparison, it is thus 
considered that an airplane, designed to fly in the speed range of interest 
to this work, is likely to ha.ve a value of Zw similar to a helicopter flying 
at the same high speed. Typical values are on the order of unity. It is 
thus concluded that the static load factor for an airplane and a helicopter 
at the same high speed is similar. 

The next approximation is a single degree of freedom approximation in 
which only vertical motion is allowed, the other two degrees (pitching and 
horizontal translation) being restrained. The body can be thought of as 
translating on a vertical rod thru its center of gravity. In this case 



Aik — A Q - O 



so, of the equations of motion, the only one that applies is 



‘H vo t£ vo 
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The aircraft response to a stead/ vertical gust is 







"floating” 



-Zwt 




normal acceleration 



This simple approximation shows that the time constant of normal 
acceleration is also dependent on Z w. In the figure are shown typical 
responses for different Zw's. 



The large Zw causes a higher peak load factor, but a faster tine 
constant in reducing the load to zero. Statically, then, a snail Zw 
is desirable, but dynamically the larger Zw is seei. to be better be- 
cause of the increased rapidity with which the gust induced accelerations 
are reduced to zero. 








The approximation obviates the fact that a lifting surface, airplane, 
or helicopter will "float" up in response to a vertical gust until the 
aircraft's vertical velocity is equal to that of the gust ar.d the angle of 
attack is thus zero. The longer it takes for this to happen, the longer 
will the aircraft experience a normal acceleration, as well as other 
forces and moments due to a non-zero angle of attack, (e.g. X K &«<). 

From the first two simple approximations it is seen that Zw, the 
lift curve slope of the aircraft, controls the load factor, and that 
automatic direct lift or thrust control could be used to either decrease 
the static load or decrease the time constant of the response, but could 
not accomplish both objectives at the same time. 

The next step is to consider two degrees of freedom, pitching and 
heaving, at constant forward speed. This is the so called short period 
approximation of the conventional, airplane, but really is an approximation 
to the motions ensuing before airspeed has a chance to change appreciably . 
The approximation is exact for the time it takes for three integrations, 
if Xw = 0. (This can be seen on the simplified analog diagram in the 
Description of Analysis Section by tracing a signal entering as w^ and 
finally causing a ). 

The equations that govern these motions are 

Cx LOc\ ~ Zv~> CsoJ 4 ~~ \J o 

and the characteristic equation is 



Zu)(S'\Ate N ) - Ht* 



l c ; 



-=■ O 



o 



It was found that Zw is likely to be of the same magnitude for the 
helicopter and the airplane, and now and must be investigated 

for contrast. The airplane will be considered first. 

For an airplane the moment derivatives depend heavily on the 
horizontal tail. The only contribution to Mq usually accounted for 
in estimates for "conventional" airplanes comes from the change in 
tail lift with tail local angle of attack as the airplane pitches. 
One way of writing the expression for Mg is: 

= it ^=1 s t P (, ..n jcn 

Se> \ i.j 



fAo 



L U. 

cl ^ 




airplane 



where Uo =U 0 ('tail) (i.e. the absence of power effects) 

and V,. = “r 1 1 Cl »>t = 

^ Sc. 



The magnitude of Mq. for the airplane thus depends on the tail size and 
the square of the tail length, and is proportional to forward speed. 
Typical values in the speed range of interest are between -1.0 and -2.0, 
as taken from the data in reference (10). 
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An expression for M ^ is easily written from a force diagram by 
summing all the moments about the center of gravity 

2 ^ ~ O ■= (xce - \o.c) Lw •+ h ; - U L t 

and taking the derivative 

= (Xce 

In the a.bsence of power effects q = q (tail) and 



kw * -Uy 



(X ci -X qQ 






For the "conventional" airplane, is thus proportional to dynamic 
pressure and depends on center of gravity position, tail length, tail 
size and fuselage characteristics. The tail term is usually the dominant 
one in the expression such that is large and negative to provide 
ample static stability. 



The "conventional" helicopter is not at all like the airplane in 
regards to and Mq . For the conventional single rotor helicopter 

with no hinge offset or horizontal tail, and neglecting the fuselage, 

Mg arises from the lag in the tilt of the resultant thrust vector as 
the shaft pitches. For these helicopters the lag of the thrust 
vector is approximately equal to the lag of the tip path plane of the 
rotor (reference 1). M ^ arises from the increased drag on the rotor 
blades (H force) as angle of attack is increased. Both of the 
derivatives depend on hub height, and are independent of center of 
gravity position. For the hub above the cenler of gravity Mg is 
negative (stable drunping term), and M ^ is positive (static instability^. 



The helicopters capable of flying in the 150 kts range, and especially 
the large, heavy ones of this study, are very unlike the "conventional” 
helicopter. These newer helicopters have main rotor flapping hinge offsets, 
horizontal tails, and large, bulky fuselages. The reason for the hinge 
offset is for greater control power, and the tail becomes necessary (and 
effective) as speed increases. This work is concerned with only these 
latter types of helicopters and all reference will be to the helicopters 
represented by the data in Appendix I. 

From the force balance diagram in Appendix (il) comes 

ZH’O = Y ^ v h) ] + U + s 

1 (1) 

+ -V tt. - JU L4 

3L 

When the small angle approximation, H = To! is made, the derivative 
taken with respect to 0, and it is noted that 

1 » , 

then 

c>H - JkT W + _ 9*. 

be Uo b 

From the relation a 1 =Ch/Ct > it can "be shown that (reference l) 



As 



C/ 



Qo/c t . 

y I t: I increases, as with heavier faster helicopters with the 



rotor plane tilted farther forward, it is seen that /b © decreases. 






It has been found characteristic of the helicopters investigated that 
this term results in a small contribution to the pitch damping at 150 kts . 
The trend of with speed for the data helicopters is shown in 

figure (3). For the 24000 lb helicopter, for instance, this term (which 
makes up all of the pitch damping in the "conventional" helicopter) 
contributes less than five percent to the total damping. 

From this information it is seen to be reasonable to neglect the tilt 
of the rotor resultant force vector in estimating the pitch damping 
at high speeds. 



The time constant of the rotor tip path plane tilting response to a 
shaft tilt rate is 



+ k. c)(X_| 



which, for the speed range and helicopters of interest can be approximated 



by 




Including these assumptions results in 
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where the ratio of the dynamic pressures in and out of the 

rotor slipstream at the tail. 
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For the data helicopters the table below shows the approximate 
breakdown for the contribution to at 150 kts. 



HELICOPTER 


©O 

(Ct/ct) 


M© 


TAIL TERi-i 


HII.'GE OFFSET TEFL'. 


24000 fwd 


3.1 


.77 


15% 


20ij 


35000 aft 


3.0 


.54 


6v?o 


35v 


42000 fwd 


2.76 


.38 


hofo 


hO% 


18000 aft 


2.97 


.53 


50% 


35% 



Typical values for the pitch damping are seen to be between -0.5 and -1.0 
at 150 kts. It is clear that the need for a tail increases with ©b /n 

A-?' 

(decreasing © ) , and the data is in agreement with the prediction of 

reference (l) which says that the tilt of the thrust vector contribution to 
the pitch damping should be zero for 



Q /(~) = 3.0 



As tail size increases the pitch damping approaches that of the airplane 
with the same tail size. 

An expression for results from taking the derivative of equation I (l'' 
with respect to - 
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Making the approximation that, for any helicopter Jk >> JU. ['Bis +- *, 

£cr/ 

using the high speed, approximate expression for /o^J derived for Zw 
earlier, and making the approximation that, in the absence of induced effects 



vJ = Uh ” 2/“* 
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then, the expression for M ^ becomes 



I (2) 
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All of the terms except Cw^ and, possibly, the second ore are 
destabilizing. As in the airplane, the tail should be large enough to 
overcome the instability of the rest of the aircraft. Unlike the airplane, 
however, there are other restrictions, such as control to trim at different 
power settings, that govern the size of the tail and have kept it relatively 
small on all present helicopters (relative to the airplane of tie sai< e gross 
weight). It can be seen from the expression above that there is little 
chance for M ^ to be negative without an extremely large tail. 






It should also be remembered that the fuselage is characteristically large 
and bulky for these helicopters, resulting in a much larger fuselage 
contribution to the static instability as forward speed increases than for 
the airplane of corresponding gross weight. The argument, from a stability 
point of view, is for a much larger tail, but as said, other factors are 
likely to keep the tail size smaller than on the corresponding airplane. 

The tilt of the resultant thrust vector term and the hinge offset term are 
both destabilizing, and both increase with a power of the forward velocity 
as can be seen from figure ( 4 ). The second term in the expression 
can be stabilizing or destabilizing, depending on the orientation of the 
resultant rotor force vector, R, with respect to the center of gravity 
in the trimmed flight condition of interest. The term is stabilizing if 
the extension of the vector, R, is behind the center of gravity, £ is 
negative, as in the figure. 
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In order for this to be the ease, the net moment in the trimmed flight 
condition contributed by the flapping of the blades about the offset 
hinges, the fuselage, and the tail must be positive. This means that the 
tail should be at negative incidence and thus carry a download to increase 
the size of this stabilizing term. A download is also desirable to in- 
crease static speed stability. Mu, but other factors such as dynamic 
stability, and control required versus trim speed and power determine the 
incidence of the tail. It should be noted that, for the "conventional" 
helicopter, neglecting the fuselage, this term is identically zero because 
the force vector, R, must pass thru the center of gravity in trim. 

For helicopters with hinge offset, tails, and large fuselages, however, this 
term is, in general, non-zero and increases with forward speed. It is also 
seen that the center of gravity position plays a smaller role in the static 
stability of the helicopter than for the a.irplane. The C.G. position, ^ 

only occurs in the <= term and its importance is obscured and diminished 
as all the other terms in equation I (2) get more important. 

The actual value of M ^ is very dependent on the particular configuration 
as each term in the expression can be of the same order of magnitude. It is 
probable, however, that will be positive or unstable for present heli- 
copters in the high speed range. The average value for the data helicopters 
at 150 kts is about +2.0. 

Returning now to the characteristic equation for the constant speed, 

A U. = 0, approximation, if "typical" poles are plotted on the complex plane, 
and the roots traced as is varied to "typical" values for the airplane 
and helicopter, the differences in characteristics become apparent. 







Helicopter 




For the helicopter, instability occurs when 



(rtoc - z„n4>) > o 



which, for "typical" values of Zw and M@ given earlier, is when 11^ 
is between +0.5 and +1.0. This is below the values for all but one of the 
data helicopters, and the effect can be clearly seen in the analog 
computer responses for the different helicopters in figure (5^. In all 
but the most stable one the motions are dominated by a rapid divergence 
in pitching - a consequence of angular acceleration induced by 
positive M^. . 

For the airplane ,the roots are of short period and well damped, 
lying in the region indicated on the figure. This means that the airplane 
will accelerate to the gust velocity very quickly, and the moments and 
forces due to angle of attack will be quickly dissipated after a 
disturbance by a steady vertical gust. 
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A further approximation to the equations of motion to adequately 
describe airplane accelerations in response to steady gusts is not needed. 
The airplane typically exhibits a phugoid, or long period motion, in 
and while = 0 (in the body axis system). The resulting accelerations 
on the airplane are small compared to the ones in the "short period" 
response. For the helicopter, however, the pitching and angle of attack 
responses are influenced by the changes in airspeed because of the nature 
of the derivative My, the forward speed static stability. For the 
"conventional" airplane, Mu = 0, but for the conventional, as well as the 
high speed large helicopter, Mu is in general non-zero. An expression 
for Mu can be written as 
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) (fuselage) and ( ) (tail) are influenced by induced 
effects or changes of the rotor slipstream with speed. It can be seen 
that the size and magnitude may depend largely on the incidence of the 
horizontal tail, if the tail is large. Positive Mu is statically stable, 
but, as will be shown later, positive Mu leads to dynamic instability . 

All of the terms in the expression vary somewhat irregularity with 
forward speed because of the induced effects of the slipstream, even at 
speeds as high as 150 kts for the data helicopters. With a large enough 
tail, at high speed, the sign of the term can be assured by adjusting 
the tail incidence, but, with a small tail, the value of Mu is very 
dependent on the particular configuration and flight conditio) (e.g. level 
or non-level flight). Of the data helicopters, the two that exhiiited the 
highest and lowest values of Mu both carried tail down loads, while t..e two 
that exhibited little change witl speed cammed a small positive or zero 
load at 150 kts. The smallest Mu's were for the smallest helicopter with 
the smallest tail. 
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An approximation that includes only this particular effect of the 
airspeed changes is presented as follows: 

It is assumed that Xu, Xw « g and Zu = 0. 

The characteristic determinant can then be written as 



SO 9 

-Mu -t-U 



= O 



and the characteristic equation, in root locus form, 

S* ( Vibes') + - O 

For a positive Mu, which is taken from the data as typical for these kinds 
of helicopters, the root loci for increasing Mu will look approximately like 
one of the following 
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From the figure it can be seen that, regardless of the au = 0 root 
locations, the helicopter will always exhibit two real convergences if 
Mu is positive. These roots can loosed be thought of as the convergences 
of angle of attack and pitching rate, the larger one being the faster angle 
of attack response. The inclusion of the other Au related derivatives is 
shown in the analog computer responses of figure (6), ’where 

a.) is for a helicopter with the "typical" derivatives (except for Zw and Mw'' 

b) is for the lowest values of all the derivatives (except Zw and Mw ) and 

c) for their highest values (except for Zw and Mw). 

It is seen that there is very little difference in the acceleration 
responses, and 0 , for the first few (=10) seconds throughout the 

range of all these derivatives. The reason for this is that the rate of 
convergence in angle of attack depends primarily on Zw, and the accelerations 
are primarily dependent on the angle of attack. For the first few seconds 
the acceleration response is almost independent of anything except Mw ar.d Zw 
and, as such, should be accurately approximated by the two degree of freedom 
model in which AIK = 0. This corresponds to the fast convergent root which 
represents the AoC motions before coupling with AW. motions. 

The result of all of the discussion above is the point that there should 
be little difference in the first few seconds of airplane and helicopter 
normal acceleration responses, regardless of the stability or instability 
of the two degree or three degree of freedom models. The static and the ore 
degree of freedom models show, with the assumption that Zw is of the same 
magnitude, that the responses 'will be identical. Inclusion of pitching as a 
degree of freedom only changes the time constant of the normal acceleration 
convergence, making it longer with positive and shorter with negative ' 

If Mo< is large and positive the convergence becomes a divergence and 
inclusion of airspeed as a degree of freedom is necessary to describe the 
motion, but when this is done (for positive Mu) it is seen that a convergerce 
reappears on the complex plane. Physical interpretation of this will be 
given in the following qualitative discussion. 
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Regardless of root locations or characteristic inodes of motion, there 
are some fundamental differences between airplane and helicopter gist 
responses that do not show up on the complex plane. The biggest, most 
obvious, and most fundamental is the fact that with different signs of 
M*, helicopters and airplanes pitch in different directions in response 
to a vertical gust. If for an aircraft were zero, its initial response 
to a vertical gust would be pure heave - no pitching at all - and a very 
slow and small response in A( A whose sign depends on Xw. (Xw can be 
positive or negative but is always small.) For the airplane with negative 
the response is a nose down pitching along with the vertical acceleration. 
Ill is causes an increase in airspeed as the important Au response (for 
Xw«g) due to the g0 term in the X force equation. The airplane then 
enters into its characteristic long period, or phugoid, oscillation. The 
opposite initial response occurs for the helicopter with positive M^. 

The helicopter pitches up, and the important response is a less of 
airspeed. In addition, if Mu is positive, as it is with all of the data 
helicopters, the AU response decreases the pitching acceleration and tends 
to cause an oscillatory mode if the effect is strong enough. The following 
diagram illustrates the differences between the two aircraft by showing what 
an observer, flying straight and level at a constant airspeed, would see 
when the two aircraft encountered gusts. 
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When the helicopter oscillatory mode is of a long enough period, the 
major part of the acceleration responses will be convergent, like the air- 
plane, and will take place in the first few seconds. A pilot would 
probably not notice the differences in the normal acceleration response 
for the two aircraft. In figure (7) are shown analog computer responses 
for the same aircraft but with 

a) positive M*, 

b) Mc< = 0, and 

c) negative M«. 

Both the constant speed approximation and the three degrees of freedom 
responses are shown, and they illustrate the observations made above. 

As a conclusion to this part of the work, it is possible to summarize 
the findings qualitatively as the differences a pilot would encounter in 
the two types of aircraft as a result of a gust disturbance. First, he would 
probably not notice the difference in normal acceleration, either the 
maximum load or the first few seconds of the dynamic response. Second, and 
most important, he would notice the different direction of rotational 
acceleration, especially if he is far ahead of the center of gravity. In 
the helicopter he will experience a local "g" response greater than that of 
the airplane because the helicopter pitches nose up. Third, he will rotice 
a smaller angular damping and a tendency for the helicopter to continue to 
rotate until the static speed stability causes restoring moments, usually 
resulting in a long period unstable oscillation. Finally, because the pilot 
is able to respond with his controls, he will notice the difference in 
control response required of him in order to dissipate the accelerations. 

In the airplane the pilot pulls his stick back to keep his nose up in 
response to a vertical gust, and in the helicopter he must push forward 
on the cyclic control to keep his nose down. 
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PART II 



ROTOR FLAPPING FEEDBACK CONTROL OF GUST INDUCED HELICOPTER DYNAMICS 



This part of the thesis pertains to the elementary stabilization of 
the heavy, high speed, single rotor helicopter which was characterized 
in the first part of the thesis. The stabilization schemes that are in- 
vestigated include a conventional rate and attitude system which is 
compared and contrasted with various types of rotor flapping, tilting 
moment, and thrust feedbacks. The impetus is to determine whether or 
not the latter are adequate to provide desirable gust response qualities 
for a helicopter. 

It is often considered desirable to reduce the amount of rotor 
flapping that occurs due to a gust disturbance. This is essentially 
an attempt at devising an "inherently stable" rotor with which to 
stabilize the helicopter. The reduced flapping is generally favorable 
to blade motion stability, but the extent to which this "improves" the 
dynamics of the helicopter is not given proper attention. Delta Three 
and Oehmichen pitch flap coupling, and variations of these, the Lockheed 
gyro system, and a hub moment sensing system are examples of these 
techniques. In most of the literature, however, only the "fixed shaft" 
dynamics of the rotor, equipped with the particular device, is studied. 
The consequence of this is that the actual effect on a helicopter's 
flying qualities is not known. In this part of the thesis a simplified 
investigation of the effects on the helicopter as would be encountered 
in actual flight are investigated. 
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Influence of Rotor Flapping on Stability Derivatives 



The fundamental problem in trying to determine ho w rotor flapping 
influences helicopter stability is to determine how much of the total 
forces and moments on the helicopter is attributable to flapping. An. 
assumption often made with the low speed conventional helicopter is 
that the rotor resultant force vector, R, is perpendicular to the 
rotor tip path plane, and that the motions of the two coincide after a 
disturbance. This assumption leads to the erroneous conclusion that 
the position, or rate of change of position, of the rotor tip path 
plane due to a disturbance or control input is more important in the 
stability and control of the helicopter than it really is. Stability 
only is considered in this work, but the same arguments as will be 
presented can be used for synthesis of helicopter control systems. 



Precisely speaking, helicopter accelerations only partially depend 
on rotor flapping, and as power , speed, and size increase, these other 
factors become of prime importance. Because of the way in which rotor 
flapping is treated in this ’work, i. e. as a dependent variable whose 
magnitude depends only on the helicopter variables, the influence of 
flapping on the stability derivatives is implicit. In order to expose 
the flapping dependent contributions to the derivatives, it is necessar,. 
to consider the flapping angles as "quasi-independent" variables and to 
take partial derivatives holding them constant. A moment derivative will 
thus look like 
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Each term in the equation has a special significance , v/iiich is listed belov;. 
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Change of moment with ( ) independent of 

rotor flapping. Fuselage, tail, tilt of the 
R vector independently of the TPP, and thrust 
changes make up this term. 



do! do, 



That component of the tilt of the R vector 
(H force) dependent on the longitudinal 
tilt of the TPP. 
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That part of the tilt of the R vector (K forced 
dependent on coning and lateral flapping angles. 

Direct hub moment due to hinge offset (or 
flapping restraint, or any other direct ca ise 
of moment due to the mechanical aspect of 
flapping as opposed to the aerodynamic forces 
caused by flapping . ) 

Zero. (all hcl'o^p^cr ’’h; nr* !rld 

const-- t \ r ’ lie f. : j-.rtinl l./lof :vc :: tnk • . x 



can be thought of as a measure of how closely the 
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rotor resultant force, R, and the perpendicular to the tip path plane, 

TPP, coincide. It is like a correlation function in that when the 
derivative is unity the R vector is perpendicular to the TPP and when 
the derivative is zero the orientation of the R vector is independent 

can be thought of as dihedral 



of the TPP. The terms and 

effect due to the coning and longitudinal blade flapping velocities,^ 



The following figure illustrates this. 
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It is seen that the longitudinal flapping velocity (lateral tilt, b, ) 
reduces the dihedral coning angle effect, thus decreasing the backward 
tilt of the R vector. It will be shown analytically in the corning 
discussion that the two effects almost cancel. 

The important distinction to be made in this analysis is the 
partial independence of a 1 and a ( , and the important result is 
"by how much?' . The angle a' can be written as 
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From the expressions above con be found the derivatives 
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They are, after much adgebra: 
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Furthermore, from the expressions derived in Part I for the moment 
balance come the following: 
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Substitution of all of the expressions derived above into equation II (l) 
yields 
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The last term in equation II (3) gives the net effect of conir.g ar.d 
lateral flapping, and can be expanded as 
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When the variable to be inserted in the ( ) is £ , the expression is 

approximately 
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but for any other variable inserted in the ( ), the expression 
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which is quite small in any condition. . The effect of lateral and 
collective flapping on the moment derivatives is thus seen to be small, 
especially when it is remembered from Part I that, for the high speed 
helicopter with a relatively large tail, most of the pitch damping 
comes from the tail. The last term in equation II (3) is considered 
zero and means that coning and lateral flapping cancel each other's 
influence on the longitudinal dynamics of the helicopter. 

The result is that the moment derivatives can be written as the 
sum of two terms, those depending on longitudinal flapping and those 
not depending on any flapping. The flapping 'dependent contribution 
can further be classified as that part due to flapping hinge offset 
and that part due to the coincidence of R vector tilt with TPP tilt. 

The latter is the term whose magnitude is often overestimated by the 
RiTPP assumption, and whose importance will be investigated in the 
coming discussion. 

It has been said that the coincidence of the R vector and the 
TPP motions is measured by 
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At very low inflows and forward velocities it is easj to see that the 
RxTPP approximation is fairly good. As speed, power, and inflow in- 
crease, however, as they do with the newer high speed helicopters, the 
sign and magnitude of the sum within the bracket is not obvious. It 
can be said, however, that in a high power climb condition X and 0o 
will both be large enough to assure that \ 



whereas in a low power descent, autorotation for instance, A and 0, 
will be small and v i , 

> 1 



In level flight it is necessary to expand the terras and make 
approximations to the expression in order to gain insight as to 
what are the important trends and contributions in terms of more 
fundamental parameters. The obvious sensitivity of the expression 
to trim conditions should instill caution when considering approximations. 
Usually good approximations a.re found to alter the magnitude and even the 
sign of the whole expression because the terms are all about the same size. 
Substituting into the bracketed part of equation II (2c) for a fl, gives 
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in which the approximation 



(j - = 1 



is made (normally < 5$, error), and into which is substituted 
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Finally, the approximation 



l !' - 



^ \ + r i/ x 



<C 1 



is made (normally < 5$ error) and the resulting expression substituted 
back into equation II (2c) to give 
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The inflow ratio, X , is 
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and is obviously negative in level flight and grows with forward speed, 

but the actual magnitude is sensitive to angle of attack and induced 

effects. Because of this the trend for a decrease in du'/X 

• 0 Q ) 



with inflow is seen, but magnitudes in level flight are still not 
obvious. The values exhibited by the data helicopters are presented 
in figure (8) and will serve to characterize the "typical" trends for 
heavy, high speed helicopters. In table II (l) are the values found 
using the exact expression given in equation II , at 150 kts. 

The sensitivity is seen to be great, for the approximations are modest. 
The approximation does, however, exhibit the right sign and trends. 
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The object of this analysis is to find what direct effects rotor 
flapping has on the stability derivatives of the helicopter, not to 
prove or disprove the validity of the RiTPP approximation. Because 
this is often a point of misunderstanding, however, it is considered 
worth while to dwell on the errors that would be incurred by the use 
of this approximation. These errors are of special concern to those 
who would devise some kind of helicopter stabilization system based 
on a reduction of rotor flapping. It will be shewn that serious 
over-estimates of "improvement" are the result of this thinking. 

When the RATPP approximation is made, equation II (l) reduces t 
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when it really should read 




where <-/ is computed with a, instead of a* . (The error incurred here 

of a 1 and a, in figure 9)* 



is indicated by the plots 



is of significance to the whole expression. 



When the terra AT V 

then the approximation can cause errors. It was seen in Part I that the 
influence of the R vector tilt in Mg is small with high inflows, or 
power settings. The approximation thus causes little change in the 
value of Mg in these conditions. As the power settings are reduced, 
of course, the R vector tilt is more important, but the approximation 
becomes better for the same reasons. The derivatives and Mu are the 
ones most affected. In these derivatives the "hT" terra is of the same 
magnitude as the other terms, and the differences in £>oJ_ and , become 

oo ~3o 

large. Figures ( 3 , 4 and 10) show the divergence between the two 
parameters, and indicate the magnitude of the errors. The static stability 
is over-estimated in both cases, the error increasing with speed. 

The approximation does not appreciably effect any of the important 
force stability derivatives. Flapping has an influence in the X force 
derivatives, Xw and Xu, because of the influence of flapping on the inplane 
H force. These derivatives are of relative unimportance to this -work and 
so the effects of flapping are considered to be only on the moment 
derivatives . 

Having now seen how flapping influences the forces ar.d moments, aid 
thus the stability derivatives, it remains to be shown by how much the 
derivatives may be changed by controlling the amount of flapping - as with 
a flapping feedback system. 



Longitudinal Flapping 



Consider first an hypothetical perfect feedback system designed to 
keep longitudinal flapping zero in some unspecified manner. If this •were 
an infinitely fast feedback loop such that a, E 0, then in the expressions 
for the moment derivatives 



The derivatives are changed by the now zero contributions from flapping 
hinge offset and the R vector tilt due to flapping. The moment derivatives 
are all reduced and the changes produce a different set of helicopter 
characteristics. It is to be noted that this is not a feedback to the 
cyclic control (which will be taken up later) but rather a computation 
of stability derivatives ignoring any contributions from longitudinal 
flapping. As such, it is now to be shown that the effective result, in 
the high inflow/speed condition, is helicopter stability and control 
characteristics similar to the same helicopter without flapping hinge 
offset. When equation II (3) is rewritten in a simple form, 




II (5) 
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it can be seen that, if the first tern in the bracket is snail with 
respect to the second, the net effect of ~ O 

is just the loss of the hinge offset term. Indeed, for the data heli- 
copters in the high speed flight condition, where is small, 

the following table shows the percentage of the first term to the 
total flapping contribution. 
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The hinge offset is included in the design of the helicopter primarily 
to increase the pilot's control power, or the amount of angular acceleration 
he can impart on the helicopter for a given stick deflection. If this 
perfect servo- system were installed, however, the pilot would lose the 
additional control power provided by the hinge offset, as well as that 
portion of the R vector tilt due to flapping. The pilot would feel 
almost as if he were flying the same helicopter but without hinge offset, 
especially for the last Wo entries in the table. This particular feed- 
back thus seems to contradict the reasoning for designing the hinge offset 
in the first place. 
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This same type of contradiction will unfortunately arise in any 
feedback system design in which there is a simple mixing of the pilot 
inputs with feedback inputs. For example, consider an aircraft auto- 

• 9 0 

pilot which can sense 8 , © , or © and feed back the signal to sum -with 
pilot inputs. If G is fed back, the pilot is limited in how much torque 
he can apply about the aircraft center of gravity, but any gust induced 
torques are also limited. Similarly, for © and © feedback, both the 
pilot and the gust disturbance are limited in how fast and how far, 
respectively, they can make the aircraft rotate. 

It has been necessary to limit the scope of this work , and thus the 
reduced control problem associated with increased stability to gusts is 
not considered. Certainly, compromises will have to be made in any 
control system such that the pilot has the final authority in the system, 
if he wants it. In this work the pilot has pushed a "hold" button and 
desires to retain an equilibrium trimmed flight condition, and never 
enters a control input. 

Along the same lines, it is very interesting to note that a negative 
a , — r- B t feedback system (opposite sign of the one just considered^ a so 
called "biased cyclic control system" , has been investigated in reference (M 
in connection with improving the control characteristics of a helicopter. 

With this negative (really positive in normal servo-convention) feedback, 
and with a suitable linkage, it was possible to reduce the control 
sensitivity (maximum pitching rate per unit stick deflection) by increasing 
the pitch damping at constant control power. It is mentioned in that 
reference that the feedback system seems to the pilot to be like adding 
hinge offset, confirming the assertion made in this work. Gust induced 
loads were not considered in that reference. It is felt that if they had 
been, the system would have been considered undesirable due to increases 
in static instability. There is little feasibility that ilat particular 
feedback system, in the high speed flight condition for hel icopte>'s with 
relatively large tails, will improve co trol charnct eristics bee- k w q 
is made up primarily by the tai] co: tribution, as was see* ii >*t I. 
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At 150 Knots the following changes in the moment derivatives of the 
data helicopters are found due to the perfect, infinitely fast feedback 
loop. 

TABLE II (3) 
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The changes in the derivatives change the helicopter characteristic 
equation, and the resulting root locations are shovm in figure (ll) for 
all of the data helicopters at 150 kts. As can be seen from the changes 
in root locations, improvement in stability is not impressive, especially 
in the aft center of gravity cases. 



The reason the improvements are not impressive is because the fuselage, 
tail, and non-flapping-dependent R vector tilt ( small) are large. 

This should be considered as typical of large, heavy, high speed helicopters, 
and will become more so as helicopter configuration approaches that of a 
compound helicopter. The truer this is, consequently, the less effective 
will be any type of helicopter stabilization system involving rotor flapping 
feedback control . 



For the one case where this feedback has succeeded ir. changing 
the sign and magnitude of M ^ enough to put the roots in an oscillatory 
region, it should be noted that the frequency of the mode is on the 
order of an airplane phugoid mode, and does not reserble the airplane 
classical "short period" mode in regard to frequency or damping ratio. 

The change in sign of the effective M ^ for both forward C.G. cases 
does change the direction of the initial pitching, however, and would 
result in a decreased "local g" response that a pilot would feel. 

The analysis presented above is for only one flight condition, 
at 150 kts, and shows that helicopter stability can be improved, however 
modestly, by reducing longitudinal flapping. It is clear that the 
improvement is a function of the relative changes in the moment derivative 
M k and M q . Both are reduced by the reduced flapping, and the new root 
locations on the complex plane depend on which one is reduced "more' . In 
other words, at each forward velocity the suppression of longitudinal 
flapping will change the stability characteristics by seme amount, ar.d 
that amount depends on the size of the flapping derivatives, ^^3(3 



and, finally, how much the rotor actually contributes to the stability 
derivative (e.g., the tail contributes the major portion of M q at high 
speeds). These will all change with speed, and the plot of figure (12) 
shows that the relative changes in the derivatives for the "average" 
helicopter are such as to lend credence to the flapping feedback as speed 
is increased. Unfortunately, it has been sham that, even at the high 
speeds where, from figure (12 ), the feedback seems most desirable, improve 
ment in stability is only modest. 



and 




the relative amount of R tilt due to 



flapping, 



Collective Flapping (Coning) 



Near the beginning of Part II it was shown that aa, and Ab, have little, 
indeed negligible, effect on the orientation of the rotor resultani force 
vector, R, except when the helicopter is pitching. In the absence of 
pitching, the ACL 0 and 4b t contributions to the h T force are almost 
equal and opposite in response to an angle of attack change. In a. static 
sense it can then be said that coning angle has no direct influence on 
the stability derivatives. It is, however, acceptable to approximate 
a proportionality of thrust and collective pitch, which in turn, in the 
absence of induced effects, is nearly proportional to the coning angle. 

It is approximated as 



It is now possible to write a feedback equation, using this simple 
approximation, that says the collective pitch of the blades is 
proportional to the coning angle 



aT = 






and, if all of the lift comes from the thrust of the rotor 




_L_ ^ ? n. _ W_ ciOu = 3— 




A 0 o a ip A" A(Xo 



If the feedback loop is infinitely fast there is some value 




of gain such that 0, and then 



cko/ 



— 0 which causes Zw = 0. 



In addition, now that the collective pitch can vary, there are changes 
in TPP tilt and R vector tilt due to the derivatives 



on the analog computer, that these changes in a, and a' are of prime 
importance when considering any kind of thrust feedback. They result 
in cross control (M $ ) coupling and can significantly alter the angular 
accelerations on the helicopter. 

For the present discussion, only the effects of reduced Zw will be 
considered. When Zw = 0 the characteristic equation of the constant 
speed approximation becomes 



where the vertical force equation now reduces to the identity, = A 0 . 
There is one unstable root for any positive Physically, this is 

the equation of motion of a weather vane. The sign and magnitude of 
M depends on the position of the pivot point. If the pivot is far 
enough back is positive, and the weather vane diverges in rotation 
when disturbed. 
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It will be shown later, in the analysis of the actual feedback systems 
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A reduction in Zw causes a decrease in damping and frequency of any 
aircraft. It also causes a reduction of initial normal acceleration. 

To the airplane , with its large negative Mj/, a reduced damping and 
frequency of the short period oscillation is usually a small M ccst ,; for 
an auto pilot designed to reduce Zw, and thus the initial lead factor. 

For a helicopter with a positive M^, however, a decreased Zw is had only 
at the expense of an increased rate of pitching divergence. With a 
positive it is mandatory to reduce the angle of attack quickly in order 
to achieve pitching stability, but this means increased load factors. 

This is clearly one of those cases in which other restrictions must be 
considered when deciding which is best. 

Two Degrees Vs. Three Degrees of Freedom 



Before going on with the analysis of any particular feedback system, 
it is of interest to see how the perfect flapping feedback systems 
considered above affect the long period oscillation, or the mode of motion 
arising when forward speed is allowed to vary. Using the approximation 
derived in Part I, the roots for the ^2000 lb helicopter have been 
computed for both the unmodified and the modified cases and are shown in 
figure (13). The root locations show a modest improvement in stability. 
The figure is not a conventional root loci, but rather the constant 
speed roots were computed for the two cases individually, and a root loci 
plotted for increasing Mu. 

It is seen from the small change in long period roots that there is 
very little effect of flapping in the long period oscillation. From this 
observation it seems reasonable to consider only the constant speed 
approximate equations of motion for a simplified analysis. It is probable 
that no improvement in helicopter stability characteristics can be made 
without improvement in the " All 20 mode", and so this approximation will be 
used in the analytical portion of the remainder of this work. ne malog 
computer results will show the complete dr-yees of* freedom, h never. 



Oehmichen and Delta Three Pitch Flap Coupling 

Flapping feedback control can be implemented on a helicopter mec baric ally 
or electronically. Linkages or sensors are used in some way to sense blade 
flapping and alter the pitch setting of the blades cyclically or collectivel; . 
There are two mechanical schemes which deserve special interest. The first 
is Oehmichen pitch flap coupling because, from an analysis of this mechanical 
device, insight into other more complex systems currently being proposed for 
helicopter stabilization can be found. The second is Delta Three pitch 
flap coupling because it has been considered as a possible means of achieving 
blade and helicopter stability for many years, but has not, to the author's 
knowledge, been investigated in connection with the dynamics of the helicopter 
as a whole. The two mechanical schemes are described below. 

With Oehmichen pitch flap coupling the cyclic variation of pitch that 
a blade "sees" (excluding any pilot inputs) is proport icnal to the flapping 
angle of a blade (in the case of a four bladed rotor, but otherwise this 
should be the angle of the TPP) 90° in front of it (in the direction of 
rotation). This is seen in the figure below for a four bladed rotor. 
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Positive Oehmichen coupling is defined here as: 
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In words, positive blade flapping resulting in a decrease in pitch of 
the blade 90° behind it (counter to the direction of rotation) is 
defined as positive. The mechanical gearing ratio is the gain H. 

The Delta Three hinge is the blade flapping hinge ar.d when the 
hinge is cocked opposite to the direction of rotation it is said that 
£>3 = ( •$. between the hinge line and the blade transverse axis). 

It is shorn positive in the figure. 



POSITIVE DELTA THREE 







With positive Delta Three, when a blade flaps up it reduces its own pitch 
angle by tan £3 . 



It is seen that, with both schemes the collective pitch is decreased 
proportional to the coning angle, but that the two devices differ in the 
cyclic component of feedback. The cyclic pitch due to the feedback, ? ^ 4 
is proportional to longitudinal flapping in the Oekrnicben system, but 
proportional to lateral flapping in the Delta Three system. 

In order to effectively compare the two systems it is convenient to 
separate the coning and cyclic components. This results in the feedback 
equations : 



Thus 
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or 
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where J = K - H, and J = K £ - tan S3 are the conventional, r echa ical 
arrangements for the two systems. 



This separation in the feedback signals cal' ’oe accomplished by suitable 
linkages or even by separate mechanical devices (reference 10, but, without 
regard for practicality, it is easier to think of the feedback as 
electrical signals originating from potentiometers mounted at the flapping 
hinges measuring blade flapping angles. In this way the signals may be 
used to drive servo-motors or hydraulic servos to vary the pitch of the 
blades at any azimuth position. The signals from different blades are 
easily combined and thus separation of the cyclic and collective corr.por.e..ts 
is possible. It is also possible with this kind of arrangement to integrate 
the feedback signals and, with suitable sensors, to create feedback signals 
proportional to the derivatives of blade flapping. 

Flapping feedback can be introduced into the equations of motion on 
the analog computer very easily when the blade flapping variables are 
treated as they are in this work. On the analog computer diagram, figure (l', 
it is seen that the flapping variables are just different combinations of the 
helicopter variables, each multiplied by a different constant. In actuality 
the different feedbacks are just unique combinations of angle of attack, 
airspeed, and pitching rate. As such, the feedback equations car. be written 
as : 




II (6) 





It should be noticed that the control variables appear on both sides 
of the equations, and there is therefore a closed loop effect which 
determines the actual "gain" of the system. This is illustrated in a 
simplified block diagram for the collective feedback. 




The portion outlined by the dotted line can be replaced by a single 
"block" which is the effective closed loop gain. One ot; er cbser'/atic 
that should be made here is that the Delta Three and Oei-nichen cyclic 
components are of opposite sign. B ^ t t 1 s increased when the rotor flaps 
back (positive flapping) with positive Oehmichcn, and , ifl y p is decreased 
when the rotor tilts to the advancing side (positive flapping'' with 
positive Delta Three. 
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The feedback equations can be rewritten in terms of effective gains 
which take into account the closed loop effects of the dependence of 
flapping on the control inputs. Derivation of these effective gains is 
shown for the Delta Three system. Equations II (6) can be rewritten as 



II (?) 



A © 0 “ 



-3 






c)oj 



cku. 






II (8) 
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Defining 



and substituting equations II (7) and II (8) into each other results ir. : 
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The Oehmichen feedback equations are analogous. T'ne-j are: 
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II (12) 
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Equations II (9) - II (12) are the complete feedback equations written 
in such o. form that they can be introduced into the equations of motion, 
equation (l), as control inputs. They are multiplied by the respective 
control derivatives (e.g. , Mg ( , etc.) and transferred to the left hard 

side of equation (l). This procedure results in what ma; be thought' of as 
modified stability derivatives. It is shown for the constant speed me: e. t 
equation modified by Oehmichen feedback. From equation (li, 

-Mw Csw) (js -tAb) 6e ~ tAvu auJq^ + Ko 0 ^6o 
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which can be modified by replacing />0 o and &B 1S b;, equations II (ll^ 
and II (12), resulting in: 



-Mw + Ke,[ e v^("’i] + Mb,^ 



Rearranging gives : 



Kvd = K\w °> l0 j 



where 



Mw - Mw + M 6l KTt;' 



Me - Mo - Me,3"[ + Mb, !<"[!' 



The procedure is identical for all of the derivatives ar.d for each 
proportional feedback system. It results in equations of motion of 
exactly the same form as without feedback, but with modified derivatives 
denoted by the (') prime. The writing of all of the modified derivatives 
is omitted here because little insight can be gained by looking at such 
complicated expressions whose terms may all be of the same size. 
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It is appropriate to analize each component of the two feedback systems 
separately before considering the complete feedback equations. There are 
three variations; lateral, longitudinal, or collective flapping may be 
used as the feedback "signal", independently of each other. 

Longitudinal (Qehmichen) Cyclic Feedback 

The feedback equation for this system is written as 



which can be substituted into the equations of motion and then, regrouping, 
results in modified stability derivatives as follows. 



positive. If Mw is positive, as for the unstabili: ed helicopter, then 
all of the primed derivatives are si' alien than tie utwdified ones. The 
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are always negative, Zg, and 





First, the "autopilot” control is now being accomplished by tilting 



the control axis in response to a change in longitudinal flapping. Ir. 
the hypothetical = 0, it should be remembered, no provision is made 
for the tilt of the control axis - the flapping is "suppressed" in sore 
unspecified manner, and any forces or moments dependent on flapping 
never even appear (see equation II (5))- In this actual feedback 

case, however, the control axis, AI'TF, is tilted as the "autopilot" response 
to flapping caused by a gust disturbance. The change of A ; IF angle of attack, 
> changes the H force, thrust, and induces flapping. All of tr.ese 

then change the moments on the helicopter. The .moment derivatives are 
changed by 



By regrouping in equation II (l4) to expose the implicit dependence on 
flapping, and substituting back into equation II (13), an expanded 
expression is found for a moment stability derivative. 
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where Mg can be expressed as: 
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It is seen that longitudinal (Oehmichen) cyclic feedback does not simply, 
relieve hub moments, but influences the moments on the helicopter due to 
each term in the above expression. It does not necessarily complete] 
relieve hub moments, whereas the hypothetical a, = 0 case, earlier 
considered, does. 



The second point of imports.r.ce is the occurrence of the cross control 
derivative, Z B) , in the expressions for the modified derivatives. This 
causes further differences between this actual feedback and the hypothetical 
Q , e 0 case. The cyclic control, , is primarily a moment producing ccr.trcl, 
but vertical forces are also produced by this control input, ar.d can change 
the characteristics of the control response if the cross control derivative, 
Z Bl , is large. An expression for Z g ( in initially level flight is 



which is a, good approximation in level flight. The derivative then becomes 



Z is hardly negligible in high speed flight where Zw is on the order of 
unity. When the approximation is substituted into the expression for Zw' , 



For all positive K the derivative is reduced, and thus t' e static non al 
loads on the helicopter will be reduced, due to cycl 'c feedback. 




but can be considerably simplified by neglecting all but the 
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In figures (lH) are shown the root locations for eac. data helicopter 
with Q,— feedback gain of unity (top part of each figure), including 
the effects of Zg )( . (The lines between roots are not conventional root 
loci.) It is seen that little difference between this and the A<?. v s O 
case shows up on the complex plane. Figures (15) and (lo) shew analog 
computer responses for the ^2000 lb helicopter with increasing gain, K, 

@ Zq ( = 0; and for varying Zg, @ constant K. Figure (l 6 ) shews responses 
for K = . 5 , and K = 1.0, with Zg ( at the va.lue for the ^2000 lb helicopter. 

The most important observation to be made concerning the analog 

• ^ 

computer responses is the change in sign and magnitude of the initial 0 
response - signifying a changed . The helicopter pitches down 
like an airplane with high cyclic feedback gain, and is thus statically 
stable 'with angle of attack. It is also possible to see that the static 
normal load factor is reduced slightly when Zp, is increased from zero. 
Finally, it is apparent that cyclic feedback tends to stabilize the 
helicopter in the longitudinal plane, but, for all of the same reasons 
above, the same feedback causes a destabilizing lateral-lcrgitudir.al 
coupling term. This is because the derivative, is positive, 

causing positive lateral flapping in response to a vertical gust. VJhei. 
positive lateral cyclic feedback is employed, /y p = K Ab| , the sign 

of the feedback is such as to increase the longitudinal flapping instead 
of reducing it> in response to a vertical gust. 



Lateral (Delta. Three) Cyclic Feedback 



The cyclic component of Delta Three pitch flap coupling -works the 
same way as the longitudinal feedback just considered, but the sign ar.d 
magnitude of the feedback "signal" are different. The sign is different 
because a rotor tilts with a positive b| due to a vertical gust, and 



This is clearly the wrong sign for longitudinal static liability. It 
is the right sign, however, for blade static stability. Right from the 
outset it is apparent that the positive Delta Three cyclic component will 
increase helicopter longitudina.1 static instability, and should be avoided 
from the point of view of this work. An investigation of negative Delta 
Three is warranted from this observation, and positive Delta Three will 
be disregarded until the collective component is included in the analysis. 

The procedure for analyzing this feedback is identical to the 
longitudinal one. The modified stability derivatives that result can be 
written as (for the constant speed approximation): 



When negative Delta Three is used, all of the derivatives are reduced. 



positive Delta Three then causes a decreased longitudinal cyclic, AR is aa. 
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It is possible to compare the two types of feedback by comparing the 
flapping derivatives. It can be shown that 
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For normal values of (11-13) it is seen that 
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and, as forward speed increases, approaches 

bod boJ 

The similarity of the two types of feedback is evident, but the reduction 
in pitch damping is smaller, and the increase in static angle of attack 
stability is greater with the negative Delta Three feedback. This is seen 
in figure (17) where analog computer traces of helicopter motions with the 
different feedbacks is shown. The positive Delta Tiwee arrangement is 
clearly destabilizing, and the negative Delta Three is clearly "better" 
than the Oehmichen feedback. In fact, the cyclic component of negative 
Delta Three feedback alone would seem to be quite a good stabilization 
device for a helicopter, on the basis of these results. 



Collective (Coning) Feedback 



The collective feedback equation can be written 

/\ ~ — 3 & do AcB^ — - ■ o 

and it should be noted that there is no distinction between positive 
Delta Three and Oehmic’nen because ABi* — 0. 

the equations of motion, as before, 

change in local linear velocity at 
the hub when pitching), the modified stability derivatives become 



When the feedback is introduced into 
and the approximation that 
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It is seen that, in addition to the anticipated result that Zw' is 
reduced, Mw is also reduced (made more statically stable) if is 
positive. An expression for Mq o can be written as 




\ x + -j ^g 1 L + eFvUjf 

bQo < 2 . ^^0 



where and * are both always positive. The sign of the 

/& 0 / 3©o 

remaining term depends on the orientation of the R vector in trim. For 
the data helicopters the second two terms are always larger than the first 
(the first term is generally negative because a negative € results in a 
stable contribution to M*), resulting in positive It can also be said 

that 3 unless the horizontal tail, fuselage, and hub contribute excessive 
nose up pitching moments in trim, must be positive and will become 
increasingly larger with forward speed. For the data helicopters Mq 
increases by an order of magnitude from oO to 150 knots. 

Positive coning angle feedback thus reduces both Zw 1 and Mw 1 for the 
helicopters represented in this study, resulting in a reduced maximum normal 
acceleration due to a gust disturbance, and an increase in angle of attack 
static stability. These changes cause changes in the helicopter dynamics as 
represented by shifts of the characteristic roots on the complex plane. 

As has been seen before, the reduced will "improve", and the decreased 
Zw* will vitiate the stability of the helicopter. The net effect will, 
of course, depend on the relative sizes of M^, Zw, 7 q. , and, especially. M q * . 
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In the following figure this point is illustrated by showing the various 
positions of the roots as Mw' and Z w are reduced, and the resulting 

positions of the roots when M@ o is large and small. 




The resulting root positions are thus very depc ident on ' 5 the laive>- V q. 

is, the greater the improvement in stability. 
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In figure (l4) (bottom part of each entry) are shorn the root 

locations for the data helicopters when J = 1. The very different 

ways in which the feedback has changed the root locations for the 

four helicopters attests to the sensitivity of the coning feedback 

to configuration. The 42000 lb helicopter exhibits the most pronounced 

change in that the helicopter is made stable with coning feedback. The 

analog computer responses of figure (l8) show the dependence of this 

helicopter on M£ q . When M^ 4 is zero, the helicopter is divergent, but 

with M set equal to its calculated value on the analog computer, the 

helicopter response is considerably imxjroved with coning angle feedback. 

This is the only helicopter of the four, however-, in which an appreciable 

improvement is shown. 

\ 

In order to determine the magnitude of the modified derivatives, 

J' as a. function of J, the mechanical gain must be known. According 
to the way this is defined, 
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the actual "closed loop" gain is different. 



The function has the form (as do all of the "closed loop" gains in 
these feedback systems) 



A 





and a sketch of this function is shown below 





The function is linear thru zero, and if B is large, an approximately 
linear region may be assumed for small A, and then A' A. 

When B is small, positive A' is "gain limited", but negative A' is 
amplified rapidly toward infinity. 
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In this case, J = A and B = 



o. o 

d6o 



5 



a,nd it is thus necessary 






investigate the nature of b<Xo/^ 
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From the expression given for a<>, 
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In high speed flight 7d Ci approaches zero (negatively) due to the reduced 
importance of induced effects. In the flight condition of interest it is 
reasonable to assume 



1 » 




b <- t J 



and a high speed estimate is 




From the sketch of the function for the "closed loop" gain it is seen 
that J' never gets-bigger than about (+) 2/3, but approaches negative 
infinity even with very small negative gain. When J' is, in fact, approach!! 
negative infinity, Zw is approaching positive infinity and the roots of the 
characteristic equation are approaching the zeroes of the equation 



(S<aX^b) - A’Zi(s-Me) = O 



when Mq = 0. This would seem to prove that negative conirg feedback is 
good, but it is only because of an omission of the effects of the coning 
feedback on tiie effective angle of attack stability, 'tv. . T ! en q )q is 
not exactly zero, Mw is approaching positive infinity and tins negative 
coning feedback is impracticable. 



In response to a vertical gust, the static response of the helicopter 



is 
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where, from equations II (9) thru II (12) and 



can be expressed as: 
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Figure (l^) has already been referred to, but new will serve as a 
summary of the possibilities of longitudinal or collective flapping 
feedback used individually. It shows that, for gains of unity, improve- 
ments in stability are not impressive, and certainly r.ot acceptable when 
the complexity of the mechanisms are considered. The analysis was r.ot 
carried beyond gains of unity because it has been shewn that there are 
blade motion stability limits at, or near, unity for both kinds of pitch 
flap coupling. References (3 and 6). These limits are the result of 
fixed shaft investigations of the dynamics "within" the rotor plane which 
have been assumed to be well damped and of high frequency in this study, 
and therefore negligible with respect to the motions of the helicopter. 

As the gains approach unity the results must be viewed with some 
apprehension as to their validity. An analysis which includes the blades 
as a degree of freedom coupled with the helicopter fuselage is needed to 
determine the exact stability limits. It is probable that these stability 
limits will not exceed those found from the simpler fixed shaft study, 
and thus gains above unity will not be investigated. 

Comparison of the Complete Feedback Equations (Delta Three and Cehmichenl 

Combinations of the collective and cyclic flapping feedbacks will now 
be analyzed with emphasis on comparing the effects of the lateral (Delta 
Three) and longitudinal (Oehmichen) cyclic components. The investigation 
is in two parts; first, a static analytical analysis which reveals hew 
a particular feedback combination changes the helicopter static stability; 
and second, analog computer responses showing the effects of the feedbacks 
on the overall helicopter dynamics. 



The signs and magnitudes of the terms .in equations II (17) determine the 
size and magnitude of the left hand side of equations II (l6) which can be 
thought of as modified static stability derivatives, Mw 1 and Zv . 

Since it is desired that static loads on a helicopter be reduced by means 
of the feedback system being employed, and since M and Z -g, , are positive, 
Z and Mg, are negative, the required signs for equations II (17) are 
and (-). It is evident upon examination of 

equations IT (17) that the two types of feedbacks will not be equally 
effective for corresponding mechanical gains. 



The collective "signals", c, 

, are different due to the 

differences between and , J" and . It has been shewn 

c>ay/ 

that kpiu /yw at high speeds, and it can be shewn from the 

definitions of j" and J^" that J" will always be larger than J^" for the 
same J. In most cases these two values will be relatively close, however, 
and figure ( 19 b) shows a typical variation of collective static feedback 
signal with gain, J, for K = -1, 0, +1. The negative values of 
correspond to reducing static normal loads and improving angle of attack 
stability (Mq o ^ ). From the figure it is seen that the cyclic 

component of feedback does not influence the collective feedback signal 
appreciably. This is because, in equations II (17b), is the 

x dominant term, even for large K or K$. Among the four' data helicopters, 
however, the 42000 lb one exhibits a marked difference in magnitude between 
J" and Jf," at high gains due to a particular configuration of small be 

(large K') and large resulting in a strong amplification 

of J” within the range of gains being considered (-1 < 3< + 1 ^ . 
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shown below. 




this helicopter high gain responses could not ever, be observed on the 
analog computer due to this high gain sensitivity. A plot of the static 
feedback signal, vs J (corresponding to figure (l9o^, for this 

"sensitive" helicopter is like the sketch below. 
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Only near the limits of useable gain, and only for this special 
configuration does this phenomenon occur, and so, in general, the static 
collective feedback signal, and thus the initial normal acceleration due 
to a gust, is relatively independent of the cyclic feedback open loop 
gain. 

In order to evaluate the effectiveness of the static cyclic feedback 
signal in changing the helicopter static stability, it is helpful to 
expand all of the terms into the more fundamental rotor aerodynamic 
derivatives. Before this is done, however, it should be noted that the 
equations II (l?a) are identical in form, of differing sign, and that 
when J = 0 and K and a,re positive, the Oehmichen feedback signal is 
stabilizing and the Delta Three is destabilizing. As J increases, however , 
the signals decrease and, for some J, change sign. After much algebra the 
expanded expressions can be written as: 
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The static cyclic feedback signals will thus change sign when 
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In the absence of induced effects, as a high speed approximation, and 
for 12 the feedbacks change sign when 
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It was shown before that J' is gain limited at about 2/3, but for J<1 (blade 
flapping stability limit), J' never will exceed about l/2, even at the 
highest speeds. With this information, it is predictable that the Delta 
Three cyclic feedback component will always be statically destabilizing, 
and that at high gains and high forward velocities the Oehmicher. cyclic 
feedback component will change from stabilizing to destabilizing, 
prediction is fairly independent of anything but airspeed. As the air- 
speed is reduced the induced effects boost the destabilizi. g limit of J', and, 
at the same time, reduce the actual value of J 1 so that at low speed there is 
no change in sign of the feedback. 

The prediction is verified in figure (l$kO where the static cyclic 
feedback signal is plotted for the Oehmichen and Delta Twee feedbacks. 

It is seen that, for the conventional mechanical arrangement of Celj ici.en 
pitch flap coupling, the static feedback signal charges to tl e "wror g" sign 
as the mechanical gain is increased to near unit;.' on the 2*1000 lb helicopter. 
The same result is found for all four data helicopters, the crossover pcint 
being close to that shown dn figure (I?'** but t 1 e slope of the culture vary I* 
with confd gurati on. The same co,-.f jpurnlioi (*h 000 1". lelicn'e - '' t a s e.-td 
a marked collective feedback sc itivit also s! c a i h d c. cl 1 ft .1' -o’, 
sensitivity, or iimrcr sod dr..-...- \1 si . -p ' . f i • e c • '■ *" 1 '. 



It looks like: 




The reason for this is the increased importance, in the cyclic feedback, 
of the increasingly larger collective feedback signals, as the mechanical 
gain of the system is increased. 

As a summary of the findings of this static analysis, it can be simply 
said that; 1) Positive collective feedback reduces static loads (normal as 
well as angular) almost independently of the sign or magnitude of the cyclic 
feedback, and Mq o plays an important part in how "good" this type of feedback 
is. 

2) Positive Delta Three cyclic component always increases the static angle 
of attack instability of the helicopter, and negative Delta Three is always 
stabilizing. 

3) Oehmichen cyclic component will improve static stability for most of the 
range of useable gain, J, but, at some high gain and high forward velocity, 
will begin to vitiate the static stability. 



r 






From these static results arises- a very interesting helicopter 
stabilization technique - negative Delta Three hinges on the main 
rotor. It has been shewn that this technique improves static stability, 
but two questions must be answered before it can be seriously considered. 

First is the question of blade motion stability (this work cannot uncover 
blade instability because of the assumptions made about the rotor plane) 
which immediately comes up because of the nature of the pitch-flap coupling 
(the pitch of a blade is increased when it flaps up) is intrinsically unstable. 
This matter has been investigated in reference (6) where it is sham that 
"The use of negative Delta Three introduces the possibility of flapping 
divergence but this is not a problem for normally required values of 
negative Delta Three” . 

The second question concerns the increase in Zw' associated with 
negative Delta Three hinges. According to the analysis of this work, 
in which only lowest order rotor plane dynamics are considered, one 
effect of negative Delta Three is to increase static normal gust loads. 

This is undesirable, but the validity of neglecting any higher order 
dynamics (gust alleviation as defined in Appendix IV) is questior able, 
especially with feedbacks such as these, and it is doubtful if the adverse 
effects of negative Delta Three are as pronounced as they seem in this 
analysis. A more complete analysis, however, is needed to verify the 
above assertion. 

The dynamic responses of one of the data helicopters with various 
feedback component combinations are presented as figures (20 - 23 ^ . 

As expected, those feedback combinations which improve static stability 
(M ^ ) are also the ones which improve the dynamic gust responses. As 
such, negative Delta Three (figures 21 and 23b) feedback is clearly the 
most promising feedback of the various ones investigated. 
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Figures (20 and 21) show that gust responses are improved for: 



OEHMICHEN 



DELTA THREE 



K = constant, J increasing 
J = constant , (K = some intermediate 



J = constant, K decreasing 
(increasing negatively) 



value before crossover to negative feedback) 

The analog computer responses are revealing and, in the case of the 
negative Delta Three feedbacks, the improvements in gust response are 
seen to be excellent. 

All of the results above are in agreement with, and may explain, the 
surprisingly successful flight test results briefly reported in reference (6). 
In that test a helicopter was fitted with negative Delta Three hinges 
(ij- -39°, K^= .81) and flown to its power limit. The investigation was 

directed at rotor blade flapping stability and amplitude , and so no data on 
the stability of the helicopter is presented in that report . 

Oehmichen and Delta Three Variations; Rate and Attitude 

Several helicopter feedback concepts that have been introduced in the 
literature can be quite adequately described by variations of the Oehmichen 
and Delta Three feedback components when only the macroscopic dynamics of 
the helicopter are being investigated. In reference (9) is proposed a thrust 
and tilting moment feedback system, in reference (13) is prese: ted the Lock- 
heed gyro system, and in reference (h) are presented various stabilization 
schemes. All of the above fit into this category. 

The proposed thrust feedback can be pretty well descriled bj the 
integral of the collective component of Oehmichen feedback, A©,> , 

where it is an adequate approximation tl at A @o A i • 
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From experience with the non- integrated feedback ( t> 0. = -j ha*'' it is 
predictable that M Q 0 plays an important part in the effectiveness of the 
integrated feedback. As can be seen from figure ( 2 4), when V.q is zero 
this feedback scheme vitiates the helicopter response, but when '■'$ is 
positive, increasing the gain increases the frequency and decreases the 
time to double amplitude of the response. This is similar to the result 
found with the non-integrated collective feedback where the static a: gle 
of attack instability was decreased due to the cross control derivative. 

The integral coning angle feedback (thrust feedback) dees not reduce 
static normal loads, and from figure (24) evidently reduces the effective 
Zw of the helicopter out of phase with the angle of attack variation 
(the integration makes the loop too "slow"), and thus only improves 
helicopter dynamic response if the cross control derivative, v @ 0 
is large and positive. 

Figure (25a) shows the results of integrating the cyclic component of 
Oebmichen feedback, B , s = kJ* 64 i dt as proposed in reference (M. 
static stability derivatives are modified with this sc. erne a. d helicopter 
dynamic gust responses are not improved. Again, the integrate or makes tie 
feedback too "slew". 



Combining the two above, or integrating the entire Oehnicher. feedback 
signal, results in the gust responses of figure (?5b). As with the 
individual component feedback signals, the loop is too "slow" ard dees rot 
improve the response of the helicopter. This particular ccmbinatio.. is 
the same as that proposed in reference ( 9)5 the Thrust and Tilting ! 
feedback . 



The one variation of tire Oebmichen feedback that deserves special 
attention at this time is the Lockheed gyro system. This system 'ns Tec 
devised in order- to relieve dynamic rotor loads a* d to stuT'lize t’<_ Heli- 
copter, as a whole, :ii the presence of giuts. Tie first fei 

treated, within the scope of tl is work, and ro the 1 Sri as d !b r Vb 

herein idealizes the evk inge c C l 1 * s r< • J ’ •• , a ? < ■ 1 „ ‘ . 1 . 

second n’llred fc>- tu.-- ■ . 



eC 1 
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In the ideal system L = 0, [f = 1 (in the nomenclature of reference 13) 
and the feedback reduces to the integral of Oehmicher. plus the integral 
of pitching rate, for the pur-poses of this work. A block diagrai. cf the 
idealized system is shown below. 
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It is not possible to perform the sere static feedback analysis as 
was done on the Oehmichen and Delta Three systems because of the i.ntegratio 
which does not "let through" any static feedback. A feedback equation can 
be written, however, which resembles the ones of the previous sections. 

From the diagram: 



At),, 
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k, 
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where = "sweep constant" = r (^~0 in reference (13) 



and C2 - linkage between gyro and s'wash 

Except for the added pitching rate term, and the lag time constant 
associated with the feedback (s in the denominator ) , this Gyro feedback 
equation resembles the Oehmichen equation. Typical values for K ^ a^e 
0 ”^l/2 and = 2/3. 



so that even with large there is a long lag time, or the feedback 
loop is "slow". Because of this long lag tine the gyro feedback dees not 
appreciably alter the motions before the angle of attack convergence, n:.d he .ce 
the normal acceleration response, is completed. Since there is no collective 
feedback there is only a small alteration in the normal force equation 
). The lagged rate feedback acts as a torsional spring, but 
slows the rotational motions only slightly. For all of these reasers 
the analog computer 4 responses (figure 25c) for the idealised Lockheed 
gyro feedback system are sorely disappointing. The responses are for 
Ci = 2 /3 and increasing K^. 

As a standard for comparison of all of these feedback techniques that 
have been looked at, it is now appropriate to show the gust response of the 
helicopter equipped with conventional SA3. Figure (26) presents the a’ alc^ 
computer responses of two of the data helicopters modified by rate mi 
attitude feedback - the last two strips of the figure are for values of 
t,, 4 ^ and with which the production models of those helicopters ar<_ 
presently equipped. 

The helicopter responses are, at first gift' ce, quit'? "good". ,t sems 
that there should bo little reason for a pilot to c>. ; ■] k- in. Ueca if 4 d. , 
however , it b-:oom 4 ‘S nectwuy to If 1 r -iv cltw 1, at 4 4 r , : 

the model used for otta' 1 t) ; r. y ■ . d, t 



It was found before that 
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on here, as this work is not designed to provide very accurate detailed 
models . 

In the figure it is seen that the acceleration responses, while hea. ilj 
damped, exhibit considerable static (time = 0 + ) values and over shoot. 
Because this feedback does not alter the static response (!■'. ), gust induced 
loads (normal a,nd torque about the center of gravity) are unchanged. This 
means that the pilot will still feel the same local "g" (i.e. the normal 
load, Zw wg, plus the upward angular acceleration times his distance from 
the C.G.) response as he did without the SAS. Even though the motions are 
then dissipated quickly, the pilot is still subject to the sane maximum 
forces on his body with this feedback. 

It is evident, from the point of view of this study, that the 
improvement in angular static stability is the prime determinant of the 
quality of a feedback technique. As such, no flapping feedback system 
with an integrator in the loop, and no body mounted SAS wbic does not 
sense angular acceleration will appreciably improve the helicopter pilot's 
opinion of the flying qualities of his craft. 
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CONCLUSIONS 



Parts I and II, being somewhat separa/te from each other in material 
presented, warrant separa/te conclusions. A summary of the findings in 
Part I is given on pages (28 and 29 ) as the physical interpretation of 
the similarities and differences found in that part. In addition: 

1) Airplane and helicopter Zw's are found to be quite nearly equ*l 
when comparing aircraft of similar high speed (150 kts) design conditions. 
This leads to the discovery tha,t the normal acceleration gust response 

of the two types of aircraft is similar. 

2) The size of Mq and the sign and magnitude of are different 
for the two types of aircraft, the helicopter being on the unstable side 
in both cases. This les,ds to the realization that the rotational motions 
are of prime importance in helicopter gust response. 

3) The influence of the fuselage, t-il, and rotor flapping hinge 
offset (or non-articulated blades) becomes more import nnt than, and in many 
cases overshadows , "conventional" helicopter stability characteristics 

in high speed flight with the "typical" helicopter characterized herein. 

This is to say that rotor flapping becomes less important in the presence 
of these other factors . 

4) Of the various stabilization techniques, involving rotor flappirg rs 
the basis for the feedback signal studied in this work, negative Delta 
(positive pitch-flap coupling) shaved the best (and only acceptable^ res Its. 
The static stability is improved considerably with thi s type of fe<vi’ do* . 

It is considered that this finding warrants further detailed stud;, in lids 
area. 



5) Integrated feedbacks, the sigials for which originate es inter 
flapping variables, are too "slow" to improve the gust rest 
of the helicopter, hecaus-' of tl.5 » t « vat r ii lb • 1 0 ■; , t. :v Is > c 
in static stabi’l it,, a 1 on ) 1.1,. . ' vl . • i 

response. 
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FIGUi!' ( 1 h ) : Coni I nut* 1 . 




FIGURE (15): Longitudinal (Oehiaichen) flapping feedback, 

B, ■ K a, , for 2^ » 0 K - .25, 0.5, .75, 1.0 (142000 lb 
helicopter) . 
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FIGURE (l6): Longitudinal (Oehmichen) flapping feedback, 

B,<i * K a,, for n) K ■ 0.5, varying cross control derivative, 
Z 8| * .o5, .25, b) Z ft) 

K » 0.5, 1.0. 
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(a) M 

FIGURu(l7): Cyclic component of Delta Three feedback for 

configuration similar to, but not exactly, the )»2000 lb hell 



copter, a) positive (conventional) , Kj ■ 0.1, b) negative 

Kg - .25, -. 75 , 1.0. 









(a) (b) 

FIGUllu (lS): Coning an^le (collective flapping) feedback, 

0 o ■ -J Op, a) )!(>„» 0, J » .25, 0.5, 1.0. b) M 0> - value of 



42000 lb helicopter, J » 0.5, 1.0. 
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FIGURE (CG). Oehmichen pitch flap coupling on U2000 lb 
helicopter , a) Effect of coning cor pone nt gain, J ■ 0, 0.5, 
1.0, t) K=*=0.5. b) Effect of cyclic conponent gain, K - .25, 



.75, 0.5. 




FIGURE (21); Delta Three pitch flap coupling. Effect of 

sign and magnitude of cyclic cosponent on configuration 

oiuilar to, but not exactly, the U2000 lb helicopter, ■ +. 25 , 

0, .25, -1.0, € J-0.5. 
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flap coupling on U2000 lb 



helicopter; K ■> J 



.25, 0.5, .75- 
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FIGURE (23): Delta Three ritch flap coupling for configuration 

siiLilar to, but not exactly, the 1*2000 lb helicopter a) Foaitive 
(conventional ) J «■ IC^ ■ -FI?* 0-5 b) ^ef atlvo J ■ » < ■ -0.5* 



. 75 , 1 . 0 . 




a) I’i£ o « vt»2ue of *42000 lb helicopter, J » .25* *75 » 1*0, 



b) ^ » 



0, j » .25 
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FIGURE (25): On the 1*2000 lb helicopter.*, a) Integrated 

cyclic component of Oehrichen feedback » < ^ * K a | ; K — • 05 > 
.25, .75. b) Intet .rated Oehnichen Feedback, B,^ ■ kj* a, dt, 
0 e « J Ja 0 dt , J ■ k ■ .25, 0.5 - , c) Idealized Lockheed 
Gyro, C 2 » 2/i. “ 0.1, O.H, 0.6. 
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FIGURC (26) : Rate and Attitude feedback; (t, s^l)6D|, * R('t, »+l)b Oj 

fa) 1*2000 lb helicopter, R - .2Gl, l) t, - 0, ti» -26, 2)1, -0, 

Iv “ 2.6, 3) t ( ».22,t x « .26, h)% - .22, t v - 2.6, b) lSOOO lb 

helicopter, R « • 55»'V * 
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- APPENDIX I 



HELICOPTER DATA 



The data contained in this appendix is taken from unpublished data on ihree 
large, high performance, single rotor helicopters. Forward and a.ft center of 
gravity positions, heavy and light gross weights, and three different configurable^ 
are represented. All of the helicopters have certain characteristics in cemor. , 
and these characteristics are taken as representative cf an; heavy, high speed, 
single rotor, helicopter. They include large flapping hinge offset (~/r = 5^), 

5 and 6 bladed main rotors, high disk loading, and relatively large horizontal 
tails . 

Some of the data did not cover the flight speed of interest ir this work, 



so it was 


interpolated 


to extend 10 knots and 15 knots in 


two cases. 


The 


following 


table shows 


some of the pertinent 


information concerning ti 


-.e heli 


copters . 
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C.G. 


MAX SPEED 


TATL 


e 


T 


WEIGHT 


LOCATION 


OF DATA (KTS) 


SIZE(s T 1,) 
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42000 


FWD. 


i4o 


1040 


5.5" 


9-3 


35000 


AFT. 


170 


1670 


5.5 ' 


8.2 


24000 


FWD. 


150 


1670 


5.5: 


8.2 


18000 


AFT. 


135 


955 


3.2b 


5.5 



The data was reduced from a space fixed axis to a stabilit. 
di mens i onali zed. At each speed the high, low and a\era ;o of eao! 
derivative was found and the inform tion is present:! in t!.c fi." 



lines thru the data points reprcsei t "t picul 



or ever a 



: I'd 



derivatives, and the trend., of \ari<'1 i* a with form \1 s] el. 



axis and 
st^bilU , 
■•or . ! e 

- 1 M; 



A I 1 
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O 17500 lbs AFT CG 
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APPENDIX III 



AXIS SYSTEMS AND AXIS CONVERSION 



There are three types of axis systems commonly used ir. aircraft 
stability analyses, each having its own peculiar advantages. The 
conversion of data from one system to another is usually required when 
data comes from different sources or when the analyst is more familiar 
with one system. The axis systems are named by the orientation of the 
X axis. First, there is the horizontal, or earth fixed axis system. 

In this system the relative wind velocities are resolved into horizontal 
and vertical components. Second, there is a wind axis system in which 
the X axis is directed into the relative wind at all times. In both 
of the axis systems above, an aircraft attitude change (e.g. angle between 
wing chord and horizon) results in a new orientation of the X axis with, 
respect to the aircraft. This changes the velocity components with 
respect to the aircraft, and moments and forces result. This results 
in non-zero stability derivatives like 

ch> X < S>H „ ^3 

2>o ' *c)o 

It is important to realize that these derivatives arise only because of 
the choice of axis system. Aerodynamic forces do not depend on pitch 
attitude of the aircraft body (e.g. wing chord) with respect to the 
horizon. 

The one axis system left is the body axis system. As the name 
implies the X axis is fixed with respect to the body (e.g. a wing chord 
or a rotor shaft) and al.l velocities are resolved parallel and 
perpendicular to the body. With such a system an attitude chrnge does 
not change the velocity components with respect to the bodjv and hence no 
aerodynamic forces arise. This system is choract er‘ d 1 ; the also uc 
of stability derivatives in 0 , but result.' ii a bit of ce fusiv. \.) • n 
trying to interpret physical.],, il r- forces n the b 1 . 
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This is because the rotating (in the case of a pitching or rolling rate') 
coordinate system introduces inertial, or centrifugal, forces in order 
to satisfy Newton's laws. Care must be taken to remember that these 
inertial forces arise only because the interpreter is in a fixed coordinate 
system and the aircraft in a rotating system. As such, the centrifugal 
forces are not really forces on the aircraft at all, but rather a way for 
the interpreter to satisfy himself that F = ma without having to switch 
to the rotating axes of the aircraft. 

A special case of the body axis system is when the X axis is initially 
aligned with the relative wind , and then fixed to the body - all at time 
zero - but then stays fixed with that orientation in the body during air- 
craft motions. This further simplifies the equations of motion by 
eliminating a wo term in the Z force equation. 

CONVERSION FROM SPACE FIXED TO STABILITY AXES 

In one way of converting data from one axis system to another the 
velocity components in the first are resolved into the second and sub- 
stituted into the equations of motion, expressed in the variables of the 
first. Regrouping results in equations of the same form but with different 
coefficients. These new coefficients are renamed the stability derivatives 
in the second axis system. This is demonstrated by the following example 
of converting from space fixed to stability axes. 



X 




v Z *t /c< 



A i • 



From the diagram of the velocity vector components comes the 
relationship 



u 







W = + U* 9 



where u, w are body axis velocities and u x , W£ are horizontal (space 
fixed) velocities. Expanding, for small perturbations, gives 



A (1) 



AU 


- aia x - 


bn. 


AO 


+ 


AW 




K, 


A0 


+ Oo AUy) 



The stability axes are defined by the initial angle of attack being 
zero so, for stability axes, 

Oo ' + Yo " 

If, in addition, initially straight level flight is considered then 

e 0 « -- o 

The velocity vector resolution, equation A (l), then becomes 

A U ^ * 

A (2) j , 

Au) = o ** + Uo/a0 

The Z force equation may be written, in a space axis, as 



A>d x 

c) IA K 






- If* i 



o 



AMI S 



Substituting the equations A (2) into the above, 






•Y AlA -}* 
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hL* Ato -v- 
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/'A © 4- 
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L 




9 



AG 

l\ x 



— G 



The first bracket term reduces to 

Gf> AG - -\J 0 t'AA6 



which is the centrifugal force that arises due to the change in axis 
systems (from stationary to rotating). The second bracketed term is 
identically zero (this is emphasized by the indication of what variables 
are held constant when those partial derivatives were taker,'). The ether 
coefficients are now renamed in the body axis as (after dividing them by mi 



A (Jt) ru C>0l A- ‘Zu + 'fcu. A K + hoc Aot - Uo AG - G 



The terms Zn and Z ^ are due to unsteady aerodynamics and are zero in Die 
work. 
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Another method that can be used to convert from one axis to 
another may be more straight forward and can provide more insight into 
the nature of the two axis systems. From the force balance diagram 
below come the expressions for body forces in terms of horizontal forces. 
(The moments are equal in both systems) 



Each body axis derivative is then found by applying the chair, rule as is 
demonstrated for Xu below. 



X = F x - F ? 0 
Z = F 2 + F k © 



M = M 




X 




•V- ^ ^ 



{&,) -F. 



?0 3u* bu 
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From equation A (l) comes 



l)lA X i 'j 


^Ux = e5 


dtA 


(boj 


<3Uy ~ ^ 






bm 


equation A (5) 




3a „ 3o ^ c 

3 (A ^ ^ U -£ 


, © = = X + C< 0 



Rewriting equation A (5) for any body 
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Finally, for stability axes 
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and 



A (6) 
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The same procedure is carried through for every derivative in the bed 
axis equations of motion. The result of all the conversion is: 



STABILITY AXIS SPACE AXIS 
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Fxu a - 


^ o [ Fxu ^ - 


Fzu x - 


i, 0 Fzu z ] 
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Mu 
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i> 0 mu 2 
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tfo Mu x 
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APPENDIX IV 



GUST ALLEVIATION 



For the simple aircraft of this work the initial (at time = 0 + ) load 
factor is instantaneous and proportional to the gust magnitude and Zw. In 
actuality, for the real helicopter or airplane, this is not the observed 
response. The differences between the actual and the model responses are 
attributed to and defined herein as gust alleviation. The definition used 
here is in contrast to another definition that would interpret this whole 
work as a study of means by which to alleviate the magnitude of a helicopter’s 
response to gusts, and thus considers artificial stabilization as gust 
alleviation. Many parameters may influence the magnitude and dynamics of 
gust alleviation. In reference (8), for example, is shown the derivation 
for a gust alleviation factor, or a corrective factor to multiply times 
the static value of Z w, which should be used to compute the initial load 
factor for a static model (e.g. rotor bolted down in a wind tunnel) due to 
a more complicated inclusion of the rotor dcwnwash. In reference (5) is 
shown the result of a detailed computer program, used to compute gust 
alleviation factors for the initial and the dynamic response by including 
unsteady aerodynamics and aero-elasticity to compute an effective value 
of Z w that varies with time and the motion of the heliccptei . A very 
detailed computer analysis such as that used in reference (5^, however, 
may or may not have included every cause of gust alleviation, and it is 
hard to gather any insight as to which are the important causes and their 
direct effects. Tne only conclusion that can really be drawn is that, for 
helicopters encountering sharp edged gusts, the effective value of Zw is 
dynamic and may vary from the static value. 

In the figure below are sham examples of what load factor responses 
could look like for the static, zero degrees of freedom model if gust 
alleviation were considered. (Note: these are not computed or predicted 

responses but are drawn for illustration only). 
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The value at time = 0 + is shown as Zw for the linear model, and KZw 
for the kind of static analysis presented in reference (8). The dotted 
line is descriptive of the results that would be obtained by extending the 
static analysis in an iterative manner - computing the thrust at each 
increment in time, finding a corresponding downwash velocity and using 
that to compute the thrust at the next point, and so on. The remaining 
line illustrates the kind of dynamic response that might be found to 
ensue with an analysis such as in reference (5). In accord with 
assumptions concerning a quasi-static analysis, and neglecting unsteady 
aerodynamics and rotor plane dynamics in this work, a constant linear 
value of Zw is used corresponding to the value to the right of the 
hatched area on the figure. 

In the static model the only aerodynamic variable that the aircraft 
responds to is AU)c^ because A© ~ Au) ~ © * 

As degrees of freedom are added the motions of the aircraft couple 
with the gust velocities to cause aerodynamic variables that are functions 
of time and the aircraft motions. As the frequency of the motions increase 
from zero (as in the static case) to the full three degrees of freedom 
motions, the errors of neglecting alleviation increase and lecre ots cured. 
The same is true if the frequency of the gust disturbances is ir creased. 
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It is helpful to think of gust alleviation as a dynamic transfer 
function to replace the Zw and It-/ potentiometers on an analog computer. 




A 



Zw 
G (jo)) 








lO 



At very low frequency the transfer function reduces to Zw, as the ' ode 
diagram shows, but at a higher frequency there may be a dynamic response. 
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Two ways of conducting an analysis, are contrasted here: A detailed 

analysis in which every conceivable cause and effect is included for 
analysis, and a simple analysis in which only the most important, or 
first order, effects are considered. This study is of the latter 
classification. Until the first order, important causes of alleviation 
can be determined and modeled simply, however, there is little justification 
for the inclusion of any alleviation in the simple study. Other investigations 
including flight testing, must be performed to determine which are the first 
order causes and what are their important effects. The scope of this paper 
does not include such investigation. 

A pertinent study would be to compare the results of simple analyses 
such as this with full scale flight test gust response data. This could 
lead to some insight as to how alleviation really influences the dynamics 
as a whole. It might then be instructive to construct different filters 
on the ana.log computer and put them between the gust inputs and the 
summing amplifiers to try to reproduce the flight test data. The 
characteristics of a particular filter that did the job might give a clue 
as to what the real causes are. 
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